Fingerprint of super-interdiffusion: anomalous intermixing in Pt/Ti 
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The ion-sputtering induced transient enhanced intermixing has been studied by molecular dynamics 
(MD) simulations in Pt/Ti and its anomalous nature has been explained as a superdiffusive transient 
enhanced interdiffusion. We find ballistic mixing and a robust mass effect in Pt/Ti. The sum of 
the square of atomic displacements ((R 2 )) asymptotically grows nonlinearily and scales as iV 2 and 
~ t 2 , where N and t are the ion- number fluence and the time of ion-sputtering, respectively. This 
anomalous behavior explains the high diffusity tail in the concentration profile obtained by Auger 
electron spectroscopy depth profiling (AES-DP) analysis in Pt/Ti/Si substrate (Pt/Ti) multilayer. 
In Ta/Ti/Pt/Si multilayer we find a linear time scaling of (-R 2 ) oc t at the Ti/Pt interface indicat- 
ing the suppression of superdiffusive features. We propose a qualitative explanation based on the 
accelerative effect of nonlinear forces provbided by the anharmonic host lattice. 
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I. INTRODUCTION 

The most well known atomic diffusion mechanisms 
(simple site exchange, vacancy, self-interstitial mediated 
or concentration dependent) are thermally activated pro- 
cesses [1]. However, there are a growing number of ev- 
idences are emerged that the anomalous broadening of 
interfaces or the high diffusity tail in the impurity con- 
centration profile could be the fingerprint of anomalously 
fast impurity diffusion [2-4]. It has been concluded by 
Abrasonis et al. that anomalously long diffusion depths 
can only be understood by the assistance of lattice accel- 
erative effects [2]. Transient enhanced diffusion has also 
been reported in post-annealed dopant implanted semi- 
conductors [5,7-11]. 

Anomalously long interdiffusion dephts has been found 
in various diffusion couples [2,3,5,7-11]. It has been re- 
ported that anomalous interdiffusion is not driven by 
bulk diffusion parameters nor by thermodynamic forces 
(such as heats of alloying) [3] or nor by heats of mix- 
ing [15,16]. Transient intermixing (IM, surface alloying) 
has also been found the most recently by molecular dy- 
namics simulations at cryogenic temperatures during the 
deposition of Pt on Al(lll) [14]. 

Computer simulations have revealed that the mass- 
anisotropy of metallic bilayers governs ion-bombardment 
induced enhancement of interdiffusion at the interface 
[16] and greatly influences surface morphology devel- 
opment [18,19] and intermixing properties during ion- 
sputtering [15,16]. In mass-anisotropic metallic bilay- 
ers, such e.g. Al/Pt, or in Ti/Pt low-energy ion- 
bombardment leads also to anomalously large broaden- 
ing at the interface [15,16,18]. The intermixing length 
(and the mixing efficiency) scales nonlinearily with the 
mass anisotropy (mass ratio) [16]. These enhanced atom- 
istic transport processes cannot be understood by the 



established mechanisms of radiation-enhanced diffusion 
(RED) [2]. This is because RED cannot cause anoma- 
lously large diffusion depths (far beyond the ion pene- 
tration depth) [2]. Moreover, the observed asymmetry of 
intermixing in Ti/Pt and in Pt/Ti can not be understood 
by the thermal spike model nor by RED [20]. Also, one 
might not explain the observed large intermixing length 
[3] and abrupt transport processes [14] by thermally ac- 
tivated diffusion. 

It has already been shown recently that the asymme- 
try of intermixing has been found in Pt/Ti and in Ti/Pt 
bilayers. However, no explanation is found for the very 
strong and much weaker interdiffusion in Pt/Ti and in 
Ti/Pt, respectively [20]. We propose in this article to 
understand the enhanced interdiffusion in Pt/Ti as a su- 
perdiffusive atomic transport process since they fulfill the 
most important condition of superdiffusion: the square 
of atomic displacements scales nonlinearily with the time 
of diffusion [32]. We present evidences in this paper for 
the fullfillment of this condition for Pt/Ti. We would 
like to show that ballistic (athermal) intermixing occurs 
in Pt/Ti in contrast to Ti/Pt in which normal radiation- 
enhanced thermally activated interdiffusion takes place. 

II. THE EXPERIMENTAL SETUP 

Although the experimental setup has been outlined in 
ref. [20], however, we briefly also give the most important 
features of the measuremnts. The experimental setup 
is the following: According to the crossectional TEM 
(XTEM) results the thickness of the layers in samples 
are: Pt 13 nm/Ti 11 nm/Si substrate (denoted through- 
out the paper as Pt/Ti), and Ta 21 nm (cap layer to pre- 
vent oxidation of Ti)/Ti 11 nm/Pt 12 nm/Si substrate 
(denoted as Ti/Pt). The XTEM images are shown in Fig 
1. For the sake of simplicity we consider our multilayer 
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samples as bilayers and we study the atomic transport 
processes at the Ti/Pt and Pt/Ti interfaces. Both sam- 
ples have been AES depth profiled by applying various 
sputtering conditions. The sample has been rotated dur- 
ing sputtering. In the following we will outline results of 
500 cV Ar + ion bombardment at an angle of incidence 
of 10° (with respect to the surface). The atomic con- 
centrations of Pt, Ti, Ta and Si were calculated by the 
relative sensitivity method taking the pure material's val- 
ues from the spectra. The oxygen atomic concentration 
has been calculated by normalizing the measured oxygen 
Auger peak-to-peak amplitude to TiC>2 [31]. The depth 
scale was determined by assuming that the sputtering 
yield (Yi) in the mixed layer is the weighted sum of the 
elemental sputtering yields (X^^^i)- 

The broadening of the interface is frequently charac- 
terizied by the depth resolution. The depth resolution 
is defined as the distance of points on the depth profile 
exhibiting 84 % and 16 % concentrations. This definiton 
has been introduced for the "normal" cases, when cither 
the ion mixing or the roughening can be described by a 
Gaussian convolution resulting in an erf fuction transi- 
tion in the depth profile. The same definition used for 
other cases as well, however. 

If the transition does differ from the erf function (e.g. 
when the mobility of one of the components of a diffusion 
couple is much higher than that of the other's) one might 
give also the distances between points of 84 % and 50 %, 
and 50 % and 16 %. The ratio of these distances gives 
us the asymmerty of intermixing (shown in Table 1). 

III. THE SETUP OF THE ATOMISTIC 
SIMULATIONS 

In order to get more insight to the mechanism of inter- 
diffusion classical molecular dynamics simulations have 
also been used to simulate the ion-solid interaction (using 
the PARC AS code [23]). Here we only shortly summarize 
the most important aspects. A variable timestep and the 
Berendsen temperature control is used to maintain the 
thermal equilibrium of the entire system. [24]. The bot- 
tom layers are held fixed in order to avoid the rotation of 
the cell. Periodic boundary conditions are imposed lat- 
crarily and a free surface is left for the ion-impacts. The 
temperature of the atoms in the outermost layers was 
softly scaled towards the desired temperature to provide 
temperature control and ensure that the pressure waves 
emanating from cascades were damped at the borders. 
The lateral sides of the cell are used as heat sink (heat 
bath) to maintain the thermal equilibrium of the entire 
system [24]. The detailed description of other technical 
aspects of the MD simulations are given in [23,24] and 
details specific to the current system in recent communi- 
cations [15,16,18,19]. 

We irradiate the bilayers Pt/Ti and Ti/Pt with 0.5 
kcV Ar+ ions repeatedly with a time interval of 10-20 



ps between each of the ion-impacts at 300 K which we 
find sufficiently long time for the termination of intcrdif- 
fusion, such as sputtering induced intermixing (ion-beam 
mixing) [15]. The initial velocity direction of the impact- 
ing ion was 10 degrees with respect to the surface of the 
crystal (grazing angle of incidence) to avoid channeling 
directions and to simulate the conditions applied during 
ion-sputtering. We randomly varied the impact position 
and the azimuth angle <f> (the direction of the ion-beam). 
In order to approach the real sputtering limit a large 
number of ion irradiation arc employed using automa- 
tized simulations conducted subsequently together with 
analyzing the history files (movie files) in each irradiation 
steps. In this article we present results up to 200 ion ir- 
radiation which we find suitable for comparing with low 
to medium fluence experiments. 200 ions are randomly 
distributed over a 20 x 20 A 2 area which corresponds to 
~ 5 x 10 15 ion/cm 2 ion fluence and the removal of few 
MLs. 

The size of the simulation cell is 110 x 110 x 90 
A including 57000 atoms (with 9 monolayers (ML) 
film/substrate). At the interface (111) of the fee crys- 
tal is parallel to (0001) of the hep crystal and the close 
packed directions are parallel. The interfacial system is 
a heterophase bicrystal and a composite object of two 
different crystals with different symmetry is created as 
follows: the hep Ti is put by hand on the (111) Pt bulk 
(and vice versa) and various structures are probed and 
arc put together randomly. Finally the one which has 
the smallest misfit strain prior to the relaxation run is 
selected. The difference between the width of the over- 
layer and the bulk does not exceed 2—3 A. The remaining 
misfit is properly minimized below ~ 6% during the re- 
laxation process so that the Ti and Pt layers keep their 
original crystal structure and we get an atomically sharp 
interface. During the relaxation (equilibration) process 
the temperature is softly scaled down to zero. Accord- 
ing to our practice we find that during the temperature 
scaling down the structure becomes sufficiently relaxed 
therefore no further check of the structure has been done. 
Then the careful heating up of the system to 300 K has 
been carried out. The systems were free from any serious 
built-in strain and the lattice mismatch is minimized to 
the lowest possible level. The film and the substrate are 
~ 20 and ~ 68 A thick, respectively. 

In order to reach the most efficient ion energy depo- 
sition at the interface, we also initialize recoils placing 
the ion above the interface by 10 A (and below the free 
surface in the 9 ML thick film) also at grazing angle of 
incidence (10° to the surface) with 500 eV ion energy. In 
this way we can concentrate directly on the intermixing 
phenomenon avoiding many other processes occur at the 
surface (surface roughening, sputter erosion, ion-induced 
surface diffusion, cluster ejection, etc.) which weaken en- 
ergy deposition at the interface. Further simplification is 
that channeling recoils are left to leave the cell and in 
the next step these energetic and sputtered particles are 
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deleted. 

We used a tight-binding many body potential, devel- 
oped by Cleri and Rosato (CR) on the basis of the sec- 
ond moment approximation to the density of states [29] , 
to describe interatomic interactions. This tight-binding 
approach is formally analogous to the embedded-atom 
method (EAM) [25] . This type of a potential gives a good 
description of lattice vacancies, including atomic migra- 
tion properties and a reasonable description of solid sur- 
faces and melting [29]. Since the present work is mostly 
associated with the elastic properties, melting behaviors, 
interface and migration energies, we believe the model 
used should be suitable for this study. The interatomic 
interactions are calculated up to the 2nd nearest neigh- 
bors and a cutoff is imposed out of this region. This 
amounts to the maximum interatomic distance of ~ 6 
A. For the crosspotential of Ti and Pt we employ an in- 
terpolation scheme [16,18,28] between the respective el- 
ements. The CR elemental potentials and the interpo- 
lation scheme for heteronuclear interactions have widely 
been used for MD simulations [16,22,26,27]. The Ti-Pt 
interatomic crosspotential of the Cleri-Rosato [29] type 
is fitted to the experimental heat of mixing of the corre- 
sponding alloy system [15,19]. The scaling factor r (the 
heteronuclear first neighbor distance) is calculated as the 
average of the elemental first neighbor distances. 

The computer animations can be seen in our web page 
[30] . Further details are given in [23] and details specific 
to the current system in recent communications [16,15]. 

IV. RESULTS 

In Fig 3 the evolution of the sum of the square 
of atomic displacements (SD) of all intermixing atoms 
(i? 2 ) = [ri{i) - Ti(t = 0)] 2 , obtained by molecu- 
lar dynamics simulations, where (rj(t) is the position 
vector of atom 'i' at time t, N is the total number of 
atoms included in the sum), can be followed as a func- 
tion of the ion flucnce. Lateral components (x, y) are ex- 
cluded from (R 2 ) and only contributions from intermix- 
ing atomic displacements perpendicular to the layers are 
included (z components). We follow during simulations 
the time evolution of (R 2 ) which reflects the atomic mi- 
gration through the interface (no other atomic transport 
processes are included). Note, that we do not calulate 
the mean square of atomic displacements (MSD) which 
is an averaged SD over the number of atoms included in 
the sum (MSD=(.R 2 }/A). MSD does not reflect the real 
physics when localized events take place, e.g. when only 
few dozens of atoms are displaced and intermixed. In 
such cases the divison by N, when N is the total number 
of the atoms in the simulation cell leads to the mean- 
ingless (R 2 )/N — > result when N — > oo, e.g. with 
the increasing number of atoms in the simulation cell. 
Also, it is hard to give the number of "active" particles 
which really take place in the transient atomic transport 



processes. Hence we prefer to use the more appropriate 
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FIG. 1. The concentration depth profile as a function of the 
removed layer thickness (nm) obtained by AES depth profiling 
analysis using ion-sputtering at 500 eV Ar + ion energy (10° 
grazing angle of incidence with respect to the surface) in P 
t/Ti (upper Fig: la) and Ti/Pt (lower Fig: lb). 

quantity SD. In Fig 3 we present (R 2 ) as a function of the 
number of ion impacts Ni (ion number fluence). (R 2 )(Ni) 
corresponds to the final value of (R 2 ) obtained during the 
Nith simulation. The final relaxed structure of the simu- 
lation of the (Ni — l)th ion-bombardment is used as the 
input structure for the A^th ion-irradiation. The asym- 
metry of mixing can clearly be seen when (R 2 )(Ni) and 
the depth profiles are compared in Ti/Pt and in Pt/Ti 
in Figs la, lb and 2. The computer animations of the 
simulations together with the plotted broadening values 
at the interface in inset Fig 3 of ref. [20] also reveal the 
stronger interdiffusion in Pt/Ti [30]. 

As it has already been shown in ref. [20] a relatively 
weak intermixing is found in Ti/Pt (a ~ 20 A) while an 
unusually high interdiffusion occurs in the Pt/Ti bilayer 
(a w 70 A). MD simulations provide ~ 8 ML (~ 20 
A) and ~ 16 ML (~ 40 A) thick interface after 200 ion 
impacts, respectively. The comparison of the measured 
a with the simulated broadening using the 84— 16 % rule 
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in both cases can only be carried out with great care [20] . 

Moreover, the applied setup of the simulation cell, in 
particular the 20 A film thickness is assumed to be appro- 
priate for simulating broadening. Our experience shows 
that the variation of the film thickness does not affect the 
final result significantly, except if ultrathin film is used 
(e.g. if less than ~ 10 A thick film). At around 5 or less 
ML thick film surface roughening could affect mixing, 
and vice versa [19]. Also, we do not carry out complete 
layer- by-layer removal as in the experiment. It turned out 
during the simulations that the ions mix the interface the 
most efficien tly when initialized ~ 10 ± 3 A above the 
interface. This value is naturally in the range of the pro- 
jected range of the ions. Hence, the most of the broaden- 
ing is coming from this regime of ion-interface distance 
. Initializing ions from the surface it takes longer and 
more ions are needed to obtain the same level of damage 
and broadening at the interface as it has been found b y 
ion-bombarding from inside the film. 

In rcf. [20] we have shown that the measured ion- 
sputtering induced broadening of a w 20 A for Ti/Pt 
is in nice agreement with the simulated value. Hence 
we expect that the most of the measured a is coming 
from intermixing and interface roughening contributes 
to a only slightly. The nice agreement could be due to 
that the saturation of intermixing (broadening) dur ing 
ion-sputtering is insensitive to the rate of mixing in the 
as received samples (a « 15 A interface width (includ- 
ing the interface roughenin g) in both samples). This is 
because a sharp interface and a weakly mixed one (be- 
fore bombardment in the as-received samples) lead to the 
same magnitude of broadening upon ion-sputtering, be- 
cause the binary systems reach the same steady state of 
saturation under the same conditions (ion energy, impact 
angele, etc.). This is rationalized by our finding that dur- 
ing simulations we start from a sharp interface (<to « 0) 
and we get a very similar magnitude of a than by AES. 

V. DISCUSSION 

A. The nonlinear scaling of (R 2 ) 

In recent papers we explained single-ion impact in- 
duced intermixing as an intcrdiffusion process governed 
by the mass-anisotropy parameter (mass ratio) in these 
bilayers [16,18]. We also present {R 2 } results for Cu/Co 
bilayer in Fig 2. Ti/Pt and Pt/Ti are mass-anisotropic, 
while Cu/Co is a mass-isotropic bilayer. Cu/Co can be 
used as a reference system, since it is composed of im- 
miscible elements, hence a weak rate of intermixing is ex- 
pected for. It has already been shown in refs. [15,16] that 
the backscattering of the hyperthermal particles (BHP) 
at the heavy interface leads to the increase in the energy 
density of the collisional cascade. We have found that the 
jumping rate of atoms through the interface is seriously 
affected by the mass-anisotropy of the interface when en- 
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FIG. 2. The simulated square of intermixing atomic dis- 
placements (R 2 ) (A 2 ) in Pt/Ti, Ti/Pt and in Cu/Co as a 
function of the ion-fluence (number of ions) obtained during 
the ion-sputtering of these bilayers at 500 eV ion energy (re- 
sults are shown up to 100 ions). The dotted lines (iso-Pt/Ti 
and iso- Ti/Pt) denote the results obtained for the artificial 
mass-isotropic Pt/Ti and Ti/Pt bilayers, respectively. The 
ions are initiated from the surface. Inset: The simulated 
broadening at the interface in A as a function of the num- 
ber of ions at 500 eV ion energy. The ions are initiated from 
10 A above the interface. The error bars denote a statistical 
uncertainty in the measure of broadening. 



crgctic atoms (hyperthermal particles) are present and 
which leads to the preferential intermixing of Pt to Ti 
[15,16]. 

According to the BHP model, in Pt/Ti, the energetic 
Ti light particles are backscattered downwards and con- 
fined below the heavy interface in the bulk which results 
in increasing energy density below the interface. In Ti/Pt 
we find the reversed case: the Ti atoms are backscat- 
tered upwards, towards the nearby free surface (the ion- 
sputtered surface is always close to the interface when 
the interface is irradiated). The injection of the Pt par- 
ticles to the Ti phase can be seen in lower Fig. 4. We 
know from ref. [15] that in Ti/Pt thermal spike occurs 
also at low ion energies. At 0.5 keV energy we also find 
thermal spike with the average lifetime of t w 5 ps aver- 
aged during few tens of consequtive events. This is due 
to the backscattering of light particles at the interface as 
it has been shown in refs. [16,15]. Hence mass anisotropy 
enhances thermal spike in both bilayers. In Pt/Ti, how- 
ever we find the further enhancement of the heat spike 
and t « 10 ps. Unfortunatelly we cannot explain further 
enhancement of intcrdiffusion in Pt/Ti by simple mass 
effect. We reach the conclusion that there must be an- 
other accelerative field which speeds up Pt particles to 
ballistic transport. This is reflected by the divergence of 
(R 2 ) from linear scaling in Fig 3. for Pt/Ti. Interdif- 
fusion takes place via ballistic jumps (ballistic mixing), 
when (R 2 ) grows asymptotically as TV 2 , where TV is the 
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number fluence (the same asymptotics holds as a func- 
tion ion-dose or ion- fluence) . This can clearly be seen 
in Fig. 3 for Pt/Ti. The horisontal axis is proportional 
to the time of ion-sputtering, hence (R 2 ) oc t 2 which is 
the time scaling of ballistic atomic transport [32]. In 
our particular case 1 ion-bombardment corresponds to 
t ~ 10 ps which we find sufficient time for the evolution 
of (R 2 )- Anyhow, above this t value the asymptotics of 
(R 2 )(t) is invariant to the choice of the elapsed time/ion- 
bombardment induced evolution of (R 2 )(t). Hence the 
transformation between ion-fluence and time scale is al- 
lowed. (R 2 ) oc t 2 and (R 2 ) oc t time scalings have been 
found even for the single-ion impacts averaged for few 
events (when (R 2 )(t) is plotted only for a single- ion im- 
pact) for Pt/Ti and Ti/Pt, respectively. 

No such ballistic behavior can be seen for Ti/Pt in Fig 

3. In Ti/Pt we find (R 2 ) oc t time scaling which is due 
to the shorter average lifetime of the collisional cascades 
[16,15]. The mean free path of the energetic particles 
are much shorter in Ti/Pt (this can be seen qualitatively 
in Fig 4. if we compare the length of the atomic tra- 
jectories for Pt between upper and lower panels of Fig 

4. ). The (R 2 ) oc t 2 scaling used to be considered as the 
signature of anomalous diffusion (superdiffusion) in the 
literature [32]. In the upper panel of Fig 2. the finger- 
print of superdiffusive feature of intermixing is detected 
by AES as a tail in the concentration profile of Pt. No 
such tail occurs for Ti/Pt in the lower panel of Fig 2. 
where the profile can be characterized by "normal" erf 
functions. Superdiffusive features, however, have never 
been reported before for intermixing, only for e.g. ran- 
dom walk of adatoms (Levy flight) on solid surfaces [32]. 
Therefore, the Pt/Ti system is highly unusual and there 
has to be more surprise in store. The anomalous impu- 
rity diffusion of N in stainless steel [2] and the observed 
large interdiffusion depths in various transition metal/ Al 
diffusion couples [3] could also be understood as super- 
interdiffusive processes. The trajectories of the reversed 
Ti recoils can be seen in upper Fig. 4 and no intermix- 
ing Ti atomic positions can be found in upper panel of 
Fig 4 (Ti/Pt). Although, Fig. 4 has no any statistical 
meaning (as it should be), however, the trajectories are 
selected among typical events hence some useful infor- 
mation can be obtained for the transport properties of 
energetic Pt atoms. In lower panel of Fig 4. we can see 
the ballistic trajectories of intermixing hyperthermal Pt 
atoms (Pt/Ti). The reversed fluxes of the Ti particles at 
the interface and the weaker intermixing of Pt atoms to 
the Ti phase result in the weaker intermixing in Ti/Pt 
than in the Pt/Ti system. This is nicely reflected in 
the (R 2 ) of Ti/Pt which remains in the range of Cu/Co. 
Hence Fig. 4 depicts us at atomistic level what we see in 
the more statistical quantity (R 2 ). In Ti/Pt we see much 
shorter inter-layer atomic trajectories while in Pt/Ti bal- 
listic trajectories of Pt atoms can be seen (going through 
the interface). 




FIG. 3. The crossectional view of a typical collisional cas- 
cade at the interface with atomic trajectories (two monolayers 
are shown at the interfaces as a crossectional slab cut in the 
middle of the simulation cell) in Ti/Pt (upper panel) and in 
Pt/Ti (lower panel). The positions of the energetic particles 
are collected up to 500 fs during a 500 eV single ion-impact 
event initialized 10 A above the interface using the same con- 
ditions as in other simulations. The vertical axis corresponds 
to the depth position given in A. The position z = is the 
depth position of the interface. 
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FIG. 4. The Cleri-Rosato potential energy (eV) as a func- 
tion of the interatomic distance (A) in Ti and in Pt averaged 
for an atom with its 12 first neighbors in the hep and fee 
lattices, respectively. 



B. The effect of mass and size anisotropy 

In order to clarify the mechanism of intermixing and 
to understand how much the interfacial anisotropy influ- 
ences IM simulations have been carried out with mass 
ratio 5 is artificially set to 5 « 1 (mass-isotropic, ions are 
initiated from the surface). We find that (R 2 ) is below 
the corresponding curve of Pt/Ti (see Fig 3, iso-Pt/Ti, 
dotted line, and animation [30]). The (R 2 ) scales nearly 
linearily as a function of the number of ions (and with 
t) for iso-Pt/Ti. The asymptotics of (R 2 ) is sensitive to 
the mass ratio. Hence we reach the conclusion that the 
mass-effect is robust, and the energetic particle scattering 
mechanism is weakened when mass-anisotropy vanishes 
and the magnitude of intermixing is also weakened. 

The same type of MD simulations has been carried 
out in Cu/Co as well, and the results are shown in Fig 3. 
The Cu/Co system is mass-isotropic, and thus we expect 
that the magnitude of (R 2 ) obtained for iso-Ti/Pt and 
iso-Pt/Ti should be similar to that of Cu/Co. It turns 
out, however, that the (R 2 ) of iso-Ti/Pt and Ti/Pt are 
higher than that of Cu/Co. The (R 2 ) of iso-Ti/Pt is 
somewhat even higher than that of Ti/Pt. Hence the 
effect of mass-isotropy on {R 2 } is opposite in Ti/Pt than 
in Pt/Ti. This might be due to the reversed fluxes of 
Ti recoils in the mass-anisotropic Ti/Pt, shown in upper 
panel of Fig 4. 

To understand the physical origin of the higher inter- 
mixing in iso-Pt/Ti than in Cu/Co, we should consider a 
more complex anisotropy of the interface then previously 
expected. We can rule out the effect of thermochemistry 
(e.g. heat of mixing AH) because it has also been shown 
recently, that AH has no apparent effect on intermix- 
ing during single-ion impact in this bilayer [15]. Also, 
thermochemistry does not explain the mixing asymme- 



try observed in Ti/Pt and in Pt/Ti since AH must be 
the same in the two systems. Moreover, simulated ion- 
sputtering with AH 3> (switching off the attractive 
term of the cross-potential and using a repulsive (Born- 
Mayer) Ti-Pt interaction potential for Pt/Ti) leads to 
intermixing instead of the expected sharp interface (see 
also animations [30]). We reach the conclusion that there 
must be a strong accelerative effect in the system which 
even counterbalance the effect of heteronuclear repulsion. 

The employed interaction potential energy functions 
are plotted in Fig (4) as a function of the interatomic 
distance for Ti and for Pt. It can clearly be seen that the 
thermal atomic sizes (the width of the potential valley at 
higher atomic kinetic energy) are different. This differ- 
ence could not be realized during the short time period of 
a collisional cascade (r < 0.5 ps). However, the increased 
lifetime of the thermal spike in these bilayers allows few 
vibrational periods and which is sufficient for the onset 
of the effect of atomic size anisotropy. The atomic radius 
(the width of the vibration around the equilibrium lattice 
site) of the Ti atoms is larger than that of in Pt. The 
atomic size anisotropy is even more robust at 1-2 eV en- 
ergies above the energy minimum (hypcrthermal atomic 
energies, see the horizontal lines in Fig 4). This explains 
the fact that the smaller and heavier Pt atoms are the 
first ballistic diffuscr during ion-beam mixing [15]. The 
atomic size difference and the mass-anisotropy together 
could explain at least partly the stronger simulated (and 
measured) intermixing in Pt/Ti and in Ti/Pt than in 
Cu/Co. In the latter mass-isotropic bilayer no serious 
atomic size difference (anisotropy) is found (see the in- 
set Fig (4)). The mixing in the iso-Pt/Ti could also be 
understood by the size anisotropy: the larger energetic 
Ti atoms are reversed at the interface and confined in 
the bulk. The smaller energetic Pt atoms, however, are 
injected to the Ti bulk during the cascade period. 



C. Nonlinear accelerative field 

Since it is evident from the analysis of the previous sub- 
sections that mass and size anisotropy do not fully acount 
for the enhancement of IM in Pt/Ti, there must be an- 
other source(s) of the divergence of (R 2 ). The nearly 
parabolic time scaling of (R 2 ) suggests ballistic IM in- 
duced by intrinsic acceleration of IM Pt particles. The 
following simple picture seems to explain the nonlinear 
time scaling of (R 2 ): Instead of the rapid solidification of 
the thermal spike intrinsic lattice effects speed up inter- 
layer transport of Pt atoms until saturation occurs in the 
Ti phase. 

Normal (thermally activated) atomic transport pro- 
cesses have already been carefully ruled out in ref. [2] 
for the anomlaous diffusion of N in nitrided austenitic 
stainless steel. There is no specific reason to expect that 
these reasonings do not hold for Pt/Ti. Namely, we rule 
out that vacancy, multivacancy or self-interstitial medi- 
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ated, stress induced or concentration dependent, thermal 
acceleration diffusion could explain the occurence of the 
asymmetry and the anomalous nature of interdiffusion in 
Pt/Ti. 

Although, the anisotropy present in the bilayer systems 
have serious effect on intermixing, the asymmetry of IM 
still remains unexplained. This seems to be puzzling at 
first sight, however, we do believe that similar nonlin- 
ear effects are responsible for the large diffusion depth of 
Pt as it has been proposed by Abrasonis et al. in ref. 
[2] . Similar nonlinear forces induced accelerative mecha- 
nism has been found the most recently for other diffusion 
couples [14]. On the basis of these findings anharmonic 
lattice vibrations (intrinsic localized modes) could pro- 
mote atomic mobility. It has already been shown that 
the localization of vibrational energy within an extended 
lattice can influence various materials properties [33]. 

Since the enchanced atomic mobility is also present in 
weakly repulsive heteronuclcar field (with repulsive Ti-Pt 
internuclcar potential) one might expect that the acceler- 
ative field is rather strong and could indeed be due to the 
intrinsic localization of anharmonic vibration modes at 
the interface region. The fundamental question remains 
to be answered yet that if nonlinear and long-range forces 
drive superdiffusion, why these nonlocal anharmonic vi- 
brational modes are active in Pt/Ti and why they are 
suppressed in Ti/Pt. Using repulsive Ti-Pt potential in 
Ti/Pt we find the further suppression of IM hence no 
accelerative long-range forces are active in this material 
which could not surpass the slow down of the particles in 
the repulsive interparticle field. 

The asymmetry of accelerative effects has also been 
found in Pt/Al couple during simulated vapor deposition 
of Pt on Al(lll). No transient intermixing is found for 
Al deposition on Pt(lll) [14]. In this case we explained 
the apparent asymmetry by atomic size mismatch [14]. 
The ultrafast diffuser specie must be the smaller one. 
Moreover the host material atoms must be not only larger 
but also anharmonic. This is indeed also the case for Ti. 
Hence anharmonic vibrational modes could be excited in 
the anharmonic host material which exhibits accelerative 
field. In the reverse case the larger particles could not 
excite the nearly harmonic host material with smaller 
particles. 

It should also be noted that the intrinsic accelerative 
effects of the anharmonic lattice modes accelerate mov- 
ing particles during the cascade and thermal spike pe- 
riod increasing the lifetime of the heat spike. Few of 
these accelerated particles could reach hyperthermal ki- 
netic energies (hot atoms) . According to our simulations 
in Pt/Ti the mean free path of these particles is in the 
range of ~ 10 A. More careful analysis certainly will 
be necessary, however, to identify the details of this spe- 
cific mechanism which could also hold for couple of other 
diffusion couples. 



VI. CONCLUSIONS 

We could reproduce by means of MD simulations 
the experimentally observed mixing asymmetry between 
Ti/Pt and Pt/Ti. We find a robust mass effect on inter- 
facial mixing in Pt/Ti, although the mass effect does not 
fully account for intermixing in Pt/Ti: the atomic size 
and mass-anisotropy are together could be responsible 
for the observed strong interfacial mixing in Pt/Ti and 
the weaker mixing in Cu/Co is explained by atomic size 
and mass-isotropy. We get a nice agreement for interface 
broadening in Ti/Pt with experiment while for Pt/Ti the 
discrepancy is not negligible. We conclude from this that 
interface roughening might has a significant contribution 
to broadening in Pt/Ti. 

We also find that the sum of the squares of atomic 
displacements ((R 2 )) through the anisotropic interface 
scales nonlinearily in Pt/Ti ((R 2 ) oc t 2 ) as a function of 
the time (and the ion-number fluence) as shown in Fig. 
3. From the upper panel of Fig. 2 we see the experimen- 
tal fingerprint of ballistic mixing (super-interdiffusion) 
as a well developed tail in the concentration profile of Pt 
with an approximate length of <~ 25 A along the depth 
direction. The long range (high diffusity) tail in the AES 
spectrum provides a compelling evidence for the superdif- 
fusive transport process of Pt atoms. In Ti/Pt a nearly 
linear scaling ((R 2 ) oc t) is found which is due to the 
shorter average lifetime of the collisional cascade period 
and which leads to shorter mean free path of the energetic 
Pt particles. The lack of a tail in the AES concentration 
profile for Ti/Pt (lower panel of Fig 2.) indicates a con- 
ccntional (thermally activated) mechanism of interdiffu- 
sion in Ti/Pt. In this system ion-sputtering increases 
broadening only slightly: in the as received sample we 
find (T ~ 15 A , and a « 20 A is measured after AES 
depth profiling. Based on atomistic simulations we find 
that the backscattering of the sizeable and light Ti par- 
ticles at the interface suppresses further intermixing. 

We conclude that the observed and simulated long 
range depth distribution of Pt atoms in the Ti phase of 
Pt/Ti cannot be understood by any established mecha- 
nisms of radiation-enhanced diffusion. In Pt/Ti, we find 
that further accelerative effects could enhance interdiffu- 
sion. A specific mechanism might come into play which 
speeds up Pt particles in the Ti bulk or at the interface. 
The divergence of (R 2 ) clearly indicates that accelera- 
tive effects are present in the lattice. At low ion energies, 
such as 0.5 keV, the damage energy caused by the slowed 
down ions cannot accelerate particles beyond the ion pen- 
etration depth. Moreover, the observed asymmetry also 
could not be understood within a simple picture of colli- 
sional cascades. The mean diffusional depth could reach 
few tens of A, which is beyond the mean free path of 
recoils at such a low ion energies. We support the expla- 
nation raised in ref. [2] that nonlinear forces could govern 
the accelereration of the particles in the bulk. However, 
it is unclear which possible scenario (breather-like soli- 
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tons, anharmonic longitudinal vibrations, etc.) explain 
the superdiffusive mechanism [2]. 
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